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Abstract: We construct random dynamics on collections of non-intersecting 
planar contours, leaving invariant the distributions of length- and area-interacting 
polygonal Markov fields with V-shaped nodes. The first of these dynamics is 
based on the dynamic construction of consistent polygonal fields, as presented in 
the original articles by Arak (1982) and Arak & Surgailis (1989, 1991), and 
it provides an easy-to-implement Metropolis-type simulation algorithm. The 
second dynamics leads to a graphical construction in the spirit of Fernandez, 
Ferrari & Garcia (1998,2002) and it yields a perfect simulation scheme in a 
finite window from the infinite-volume limit. This algorithm seems difficult to 
implement, yet its value lies in that it allows for theoretical analysis of thermo- 
dynamic limit behaviour of length-interacting polygonal fields. The results thus 
obtained include the uniqueness and exponential a-mixing of the thermodynamic 
limit of such fields in the low temperature region, in the class of infinite-volume 
Gibbs measures without infinite contours. Outside this class we conjecture the 
existence of an infinite number of extreme phases breaking both the translational 
and rotational symmetries. 
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1 Introduction 



An example of a planar Markov field with polygonal realisations was first intro- 
duced in Arak (1982). The original Arak process in a bounded open convex set 
D is constructed as briefly sketched below. We define the family Fd of admis- 
sible polygonal configurations on D by taking all the finite planar graphs 7 in 
D U dD, with straight-line segments as edges, such that 

(PI) the edges of 7 do not intersect, 

(P2) all the interior vertices of 7 (lying in D) are of degree 2, 

(P3) all the boundary vertices of 7 (lying in dD) are of degree 1, 

(P4) no two edges of 7 are colinear. 

In other words, 7 consists of a finite number of disjoint polygons, possibly 
nested and chopped off by the boundary. Further, for a finite collection (Z) = 

('i)"=i straight lines intersecting D we write To{l) for the family of admissible 
configurations 7 with the additional properties that 7 C U and j OU is a 

single interval of a strictly positive length for each li,i = 1, ...,n, possibly with 
some isolated points added. Let be the restriction to D of a homogeneous 
Poisson line process A with intensity measure given by the standard isometry- 
invariant Lebesgue measure /i on the space of straight lines in K^. One possible 
construction of /z goes by identifying a straight line I with the pair (0, p) G 
[0,7r) X R, where (psin((/)), /9cos(^)) is the vector orthogonal to I, and joining 
it to the origin, and then by endowing the parameter space [0, tt) x M with the 
usual Lebesgue measure. With the above notation, the polygonal Arak process 
Ad on D arises as the Gibbsian modification of the process induced on by 
Ad, with the Hamiltonian given by the double total edge length, that is to say 

p.. ^r)- ^^7€r.(AB)nG6^P(-21cngth(7)) 
n^D - EE,er.(A.) exp(-21ength(7)) 

for all G C Yd Borel measurable, say with respect to the usual Hausdorff dis- 
tance topology, sec Section 4 in Arak & Surgailis (1989). The Arak process has 
a number of remarkable properties. It is exactly solvable (an explicit formula 
for the partition function is available), consistent {Ad coincides in distribu- 
tion with the restriction of Ac to £> for C 2 D) and enjoys a two-dimensional 
Markov property stating that the conditional behaviour of the process in an 
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open bounded domain depends on the exterior configuration only through arbi- 
trarily close neighbourhoods of the boundary, see ibidem. These nice features 
are shared by a much broader class of processes, so-called consistent polygo- 
nal Markov fields, introduced and investigated in detail in Arak & Surgailis 
(1989, 1991). Arak, Clifford & Surgailis (1993) introduce an alternative point- 
rather than line-based representation of these models. Our description below 
specialises for the standard Arak process Ad- For a given point configuration 
X = {xi, ...,Xn} C Z) U dD denote by Tnix) the family of admissible config- 
urations 7 whose vertex set coincides with x. Write Hd for the Poisson point 
process in _D U dD with the intensity measure given by the area clement on D 
and by the length element on dD. By Theorem 1 ibidem (see also (2.6) there) 
the Arak process Ad coincides with the Gibbsian modification of the process 
on Td induced hy Hd with the Hamiltonian 

$(7) := 2 length(7) ^ loglength(e) - ^ log|sin0j,|, (2) 

where E(j) and V{"/) are, respectively, the edge and vertex sets of 7 while (p^ 
stands for the angle between the edges meeting in x if x € D and for the angle 
between the edge and the tangent to dD at x if a; e dD. This means that 

lEyi^cr m ^nr-exp(-$(7)) 
F{Ad eG)= ^gr.(nz.)ng "^' ^'^ (3) 
IEE^er^(n^)exp(-$(7)) 

for all Borel G C Td- The third equivalent description of polygonal Markov 
fields is available in terms of equilibrium evolution of one-dimensional particle 
systems, tracing the polygonal realisations of the process in two-dimensional 
time-space. This description, usually referred to as the dynamic representation 
and introduced already in the original Arak work (1982), turned out to be 
very useful in establishing the essential properties of the models. Below, we 
discuss the dynamic representation for the Arak process, see Section 4 in Arak 
& Surgailis (1989). We interpret the open convex domain D as a set of time- 
space points {t,y) £ D, with t referred to as the time coordinate and with y 
standing for the spatial coordinate of a particle at the time t. In this language, 
a straight line segment in D stands for a piece of the time-space trajectory 
of a freely moving particle. For a straight line I non-parallel to the time axis 
and crossing the domain D we define in the obvious way its entry point to 
D, in{l, D) 6 dD and its exit point out(/, D) e dD. 

We choose the time-space birth coordinates for the new particles according 
to a homogeneous intensity tt Poisson point process in D (interior birth sites) 
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superposed with a Poisson point process on the boundary (boundary birth sites) 
with the intensity measure 

k{B) = Ecard{/ e A, in(/, D) e B}, B C dD. (4) 

Each interior birth site emits two particles, moving with initial velocities v' and 
v" chosen according to the joint distribution 

e{dv',dv") ■.= t:-^\v' -v"\{l + v'^)-'"\l + v"^)-'"^dv'dv" . (5) 

This can be shown to be equivalent to choosing the directions of the straight 
lines representing the space-time trajectories of the emitted particles according 
to the distribution of the typical angle between two lines of A, see Sections 3 and 
4 in Arak & Surgailis (1989) and the references therein. Each boundary birth 
site X € dD yields one particle with initial speed v determined according to the 
distribution 9x{dv) identified by requiring that the direction of the line entering 
D bX, X and representing the time-space trajectory of the emitted particle be 
chosen according to the distribution of a straight line / G A conditioned on the 
event {x = in(Z, -D)}. 

All the particles evolve independently in time according to the following 
rules. 

(El) Between the critical moments listed below each particle moves freely with 
constant velocity so that dy ~ vdt, 

(E2) When a particle touches the boundary dD, it dies, 

(E3) In case of a collision of two particles (equal spatial coordinates y at some 
moment t with {t, y) & D), both of them die, 

(E4) The time evolution of the velocity Vt of an individual particle is given by 
a pure-jump Markov process so that 

^{vt+dt G du\vt = v) = q{v, du)dt 

for the transition kernel 

q{v, du) := \u - -^1(1+ v?)~^''^ dudt. 

It has been proven (see e.g. Lemma 4.1 in Arak & Surgailis (1989)) that with the 
above construction of the interacting particle system, the time-space trajectories 
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traced by the evolving particles coincide in distribution with the Arak process 
Ad- Moreover, a much broader class of consistent polygonal Markov fields admit 
analogous dynamic representations, possibly enhanced to allow vertices of higher 
degrees (3 and 4), see ibidem. The question of characterising the class of all 
polygonal Markov fields admitting dynamic representation is far from being 
trivial and a conjectured description of this class has been provided in Arak, 
Clifford & Surgailis (1993). 

The above dynamic construction of the Arak process makes it very suitable 
for simulation. However, in the present paper we focus our interest on the 
family of processes .4^'''', a,/? G M, arising as the Ising-like length- and area- 
interacting Gibbsian modifications of Ad- To this end we colour the original 
Arak process Ad as follows. Requiring that the polygonal contours of Ad stand 
for interfaces between black- and white-coloured regions in D leaves us almost 
surely with two possible ways of colouring D in black and white, arising from 
each other by a simple colour flip. We choose one of these colourings at random, 
with probability 1/2, thus obtaining a coloured version of ^d, denoted in the 
sequel Ad- The family of all admissible coloured polygonal configurations in 
D, carrying information not only about the planar contours it consists of, but 
also about the associated colouring, will be denoted by I'd. With this notation 
and terminology we define the (coloured) processes by 

dC{AD) Eexp(-wL-^l(i,))'^ ^ ^ 

with C{-) denoting the law of the argument random object and with 

^^'"^(7) := aA(black[7]) + /31ength(7), (7) 

where black[7] is the black-coloured region in D for 7 while A{-) stands for the 
area measure. We also write for the contour ensemble of with 

the colours 'forgotten' and, likewise, 7 for the colourless version of 7 S Td- 
Note that using the symmetry between black and white and possibly flipping 
the colours, whenever convenient we may assume without loss of generality that 
a > (and we do so in the proof of Theorem 3 below). 

Observe that the modifications of the type (6) fall into the general setting 
considered by Arak & Surgailis (1989) only for /? > 0, see Corollary 4.1 there. 
However, we find it natural to admit also negative /3's since there is no obvious 
infinite temperature non-interacting field available as the reference object for 
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polygonal Markov fields. Consequently, in the sequel we will abuse the language 
by referring to large positive values of (3 as to the low temperature region, and 
to small, possibly negative (3's as to the high temperature regime. For f3 < 
one has to check that the partition function Eexp ^— 7Y^''^'(.Ad)^ is finite. In 
Corollary 2 wc show that this is indeed the case and, consequently, the definition 
(6) is correct for all /3 G IR.. Clearly, there are no such problems for a, since the 
overall black or white area is deterministically bounded by A{D). It should be 
emphasised though that at present we are able to establish the existence of the 
thermodynamic limit only for (3 > 0, sec Theorem 3. 

Models of the type (6) have recently found interest in the physical literature, 
see Nicholls (2001). In particular, it has been argued that they exhibit a phase 
transition similar to that of the planar Ising model, with the low temperature 
phase admitting only finite contour nesting (as rigorously shown in Nicholls 
(2001)), and with the high temperature phase conjectured (not yet proven) to 
exhibit infinite contour nesting. 

Below, we shall also consider versions of the above models with empty bound- 
ary conditions, arising by conditioning the original model on the event of there 
being no vertices on the boundary, so that 

c := c 1 A^S'^^ ndD = <Dy (8) 

In particular. 

Likewise, we shall consider versions of these models with black (or white) bound- 
ary conditions given by 

^ (-^Di'bLcwhite)) ■■= ^ (A""' I -4^'''' ndD = ?,, dD is black (white) ) 



with 



A — i[°'Ol 

•^£)|black(white) •— -^D |black( white)" 



As a direct conclusion from (6) we get 

d£(i^"fl) exp(-W^"''^l(7)) 

dCiAn\k6) Eexp(-W^''^l(^D|bd)) (10) 

for bd e {0, black, white}. Observe that, unlike the unconditioned finite- volume 
fields A^l^'^\ a ^ 0, the conditioned fields with monochromatic boundary con- 
ditions are well defined also for non-convex bounded open D with piecewise 
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smooth boundary. Indeed, take any bounded open convex set D' containing D 
and set -4^|'^j'j, bd G {black, white}, to coincide with J^^f"^ conditioned on the 
event that no edge hits dD and that the colour on dD agrees with that specified 
by bd. The Markov property of polygonal fields (see Arak & Surgailis (1989)) 
implies that this construction does not depend on the choice of D' . Note that 
this argument does not apply for the empty boundary condition bd = 0, unless 
a = 0. 

The purpose of this paper is to construct for a,/? S R a family of ran- 
dom dynamics on f ^ which leave the distribution of invariant. This 
yields simulating algorithms for both of the Metropolis type and of per- 
fect type in the spirit of Fernandez, Ferrari & Garcia (1998,2002). While the 
Metropolis algorithm is given for all a, /3 e M and can be readily implemented 
(which is a subject of our work in progress), the perfect scheme is restricted 
to a = and seems more difficult to implement, yet its value lies mainly in 
that it provides important theoretical information about the thermodynamic 
limit behaviour of yll^'^'l in the low temperature region (for large fi) and in that 
it can be used to simulate in finite windows directly from the thermodynamic 
limit. The finite volume dynamics are discussed in the next Section 2. In Sec- 
tion 3 we discuss infinite-volume thermodynamic limits of polygonal fields and 
establish their existence. For q = and (3 large enough one of our dynamics, 
constructed in Subsection 2.2 below, admits an infinite-volume extension and, 
as mentioned above, it yields a perfect simulation scheme which enables us to 
show in Section 4 that for there exists exactly one thermodynamic limit 
without infinite chains, as made specific below, and that this limit is isomctry 
invariant as well as exponentially a-mixing. In particular, it follows that in 
the class of infinite-volume measures without infinite chains there exist exactly 
two extremal infinite- volume Gibbs measures for the black-dominated 
and white-dominated phase, corresponding to the same contour distribution. In 
this context it should be noted that this simple picture does not seem to extend 
to the whole simplex of infinite-volume Gibbs measures for ^[O''^!: we conjec- 
ture the existence and sketch, in Section 3 below, a tentative construction of an 
infinite number of infinite-volume states admitting infinite chains and breaking 
both the translational and rotational symmetry. 

As already mentioned above, the implementation of the algorithms described 
in this paper is a subject of our current work in progress. It should be empha- 
sised that an algorithm for simulating polygonal Markov fields, very different 



7 



than ours, has aheady been given in the hterature by Chfford & Nichohs (1994). 

2 Finite volume dynamics 

Below we construct two famiUcs of random dynamics which leave invariant the 
laws of the Gibbs-modified polygonal random fields ^^''^^ in a bounded open 
convex domain Z) C M^. First of these dynamics, leading to a practically feasible 
and easy to implement Metropolis-type simulation algorithm, is based on the 
dynamic representation of the Arak process. The second one relies mainly on 
the point- and line-based representation of general polygonal Markov fields and, 
after some additional work, leads to a graphical construction and a perfect 
algorithm discussed in Section 4. We postpone the proof of the finiteness of the 
partition function in (6) to Corollary 2 below. 

2.1 Disagreement loop birth and death dynamics 

An important concept below will be that of a disagreement loop, borrowed from 
Schreiber (2004), Section 2.2. This arises from the dynamic construction of the 
Arak process as provided by the evolution rules (El-4) with the corresponding 
birth rules, see (4) and (5). 

Suppose that we observe a particular realisation 7 G of the colourless 
basic Arak process Ad and that we modify the configuration by adding an 
extra birth site xq to the existing collection of birth sites for 7, while keeping 
the evolution rules (El-4) for all the particles, including the the two newly 
added ones if € f and the single newly added one if xq € dD. Denote the 
resulting new random (colourless) polygonal configuration by 7 © xg. A simple 
yet crucial observation is that for xq £ D the symmetric difference 7A[7®xo] is 
almost surely a single loop (a closed polygonal curve), possibly self- intersecting 
and possibly chopped off by the boundary. Indeed, this is seen as follows. The 
leftmost point of the loop 7 A [7® a;o] is of course xq ■ Each of the two new particles 
Pi,P2 emitted from xq move independently, according to (El — 4), each giving 
rise to a disagreement path. The initial segments of such a disagreement path 
correspond to the movement of a particle, say pi, before its annihilation in the 
first collision. If this is a collision with the boundary, the disagreement path 
gets chopped off and terminates there. If this is a collision with a segment 
of the original configuration 7 corresponding to a certain old particle p.-j, the 
new particle pi dies but the disagreement path continues along the part of the 
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trajectory of ps which is contained in 7 but not in 7 © xq. At some further 
moment ps dies itself in 7, touching the boundary or kiUing another particle 
P4 in 7. In the second case, however, this collision only happens for 7 and 
not for 7 ® aJo so the particle P4 survives (for some time) in 7 © yielding a 
further connected portion of the disagreement path for pi, which is contained 
in 7 © a;o but not in 7 etc. A recursive continuation of this construction shows 
that the disagreement path initiated by pi consists alternately of connected 
polygonal subpaths contained in [7 © a;o] \ 7 (call these positive parts) and in 
7 \ [7 ffi xo] (call these negative parts). Note that this disagreement path is self- 
avoiding and, in fact, it can be represented as the graph of some piecewise linear 
function 1 1-^ y{t)- Clearly, the same applies for the disagreement path initiated 
by p2- An important observation is that whenever two positive or two negative 
segments of the two disagreement paths hit each other, both disagreement paths 
die at this point and the disagreement loop closes (as opposed to intersections 
of segments of distinct signs which do not have this effect). Obviously, if the 
disagreement loop does not close in the above way, it gets eventually chopped 
off by the boundary. We shall write A®[a;o;7] = l^ll ® ^o] to denote the 
(random) disagreement loop constructed above. It remains to consider the case 
Xo S dD, which is much simpler because there is only one particle emitted and 
so A®[xo;7] = 7^ [7 ® ^n] is a single self-avoiding polygonal path eventually 
chopped off by the boundary. We abuse the language calling such A®[xo;j] a 
(degenerate) disagreement loop as well. 

Likewise, a disagreement loop arises if we remove one birth site xq from the 
collection of birth sites of an admissible polygonal configuration 7 G F/j, while 
keeping the evolution rules for all the remaining particles. We write 7 xq for 
the configuration obtained from 7 by removing xq from the list of the birth sites, 
while the resulting random disagreement loop is denoted by A®[a;o;7] so that 
A©[a;o;7] =7A[7ea;o]. 

With the above terminology we are in a position to describe a random dy- 
namics on the coloured configuration space Td, which leaves invariant the law 
of the basic Arak process Ad- Particular care is needed, however, to distinguish 
between the notion of time considered in the dynamic representation of the 
Arak process as well as throughout the construction of the disagreement loops 
above, and the notion of time to be introduced for the random dynamics on Tjj 
constructed below. To make this distinction clear we shall refer to the former as 
to the representation time (r-time for short) and shall keep for it the notation 
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t, while the latter will be called the simulation time (s-time for short) and will 
be consequently denoted by s in the sequel. 

Consider the following pure jump birth and death type Markovian dynamics 
on tn- 

(DL:birth) With intensity [wdx + K{dx)]ds set js+ds := 7s ® a; for k as in (4), 
then construct js+ds by randomly choosing, with probability 1/2, either 
of the two possible colourings for "fs+ds, 

(DL:death) For each birth site x in 7^ with intensity 1 set js+ds '■= 7s O x, 
then construct js+ds by randomly choosing, with probability 1/2, either 
of the two possible colourings for 7s+ds • 

If none of the above updates occurs we keep 7s+ds = 7s- It is convenient to 
perceive the above dynamics in terms of generating random disagreement loops 
A and setting 7s+ds := 7s AA, with the loops of the type A®[-, •] corresponding 
to the rule (DL:birth) and A©[-, •] to the rule (DL:death). 

As an direct consequence of the dynamic representation of the Arak process 
we obtain 

Proposition 1 The distribution of the Arak process Ad is the unique invariant 
law of the dynamics given by (DL:birth) and (DL:death). The resulting 
stationary process is reversible. Moreover, for any initial distribution of 70 the 
laws of the random polygonal fields 7^ converge in variational distance to the 
law of Ad as s ^ 00. 

The uniqueness and convergence statements in the above proposition require a 
short justification. They both follow by the observation that, in finite volume, 
regardless of the initial state, the process % spends a non-null fraction of time 
in the state 'black' (no contours, the whole domain D coloured black). Indeed, 
this observation allows us to conclude the required uniqueness and convergence 
by a standard coupling argument. 

Below, we show that the laws of the Gibbs-modified polygonal fields A^^'^^ 
arise as the unique invariant distributions for appropriate modifications of the 
reference dynamics (DL:birth), (DL:death). The main change is that the 
birth and death updates are no more performed unconditionally, they pass an 
acceptance test instead and are accepted with certain state-dependent probabili- 
ties, upon failure of the acceptance test the update is discarded. For a > 0, 6 > 
and a + a>0,f3 + b>0 consider the following dynamics 
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(DL : birth[a, /3;a, b]) With intensity [ndx + K{dx)]ds do 

• put 5 := 7s © X, 

• construct S by randomly choosing, with probabihty 1/2, either of the 

two possible colourings for S, 

• accept 6 with probability 

exp (^-aA (black[,5] \ black[7s]) - /31ength(<5 \ 7s)) 

exp (^-aA (^black[i5]Ablack[7s]) - folength(5A7s)) , 

• if accepted, set 7s+ds := S, otherwise keep -js+ds '■= 7s- 
(DL : death[a, /3; a, b]) For each birth site x in 7s with intensity 1 do 

• put (5 := 7s X, 

• construct S by randomly choosing, with probability 1/2, either of the 
two possible colourings for S, 

• accept 6 with probability 

exp (^-aA ^black[5] \ black[7s]) - /?length(5 \ 7s)) 

exp (^-aA (^black[^]Ablack[7s]) - 61ength((5A7s)) , 

• if accepted, set 7s+ds '■= otherwise keep 7s+ds '■= 7s- 

In analogy with its original reference form (DL:birth), (DL:death), the above 
dynamics should be thought of as generating random disagreement loops A and 
setting 7s+ds ■— 7AA provided A passes the acceptance test. It should be 
emphasised that the random disagreement loops above are generated according 
to the dynamic representation of the original Arak process Ad- The following 
theorem justifies the above construction. 

Theorem 1 For each a > 0, 6 > and a + a > 0,(3 + b > the law of the 

Gibbs-modified Arak process i^''^! IS the unique invariant distribution of the dy- 
namics (DL : birth[Q!, /?; a, b]), (DL : death[Q!, /?; a, b]). The resulting station- 
ary process is reversible. For any initial distribution of% the laws of the random 
polygonal fields converge in variational distance to the law ofA^^''^^ as s ^ 00. 
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Theorem 1 can be easily concluded from Proposition 1 by a straightforward 
check of the detailed balance conditions. We chose, however, to provide below 
a geometric proof of this result for the case a,(3 >0, revealing, in our opinion, 
the geometric intuition underlying the dynamics (a similar proof can be pro- 
vided for a < or /3 < as well). Note that the reason for introducing the 
additional parameters a and b with the possibility that a>0, 6>0, Q;-|-a>0 
and p + b> was to gain direct control over the diameter of the region affected 
by a single update, which decays exponentially in the current dynamics. The 
control of the diameter of the affected region is a condition sine qua non for 
possible infinite volume extensions of the (DL : ...[a, /?; a, b]) dynamics, which 
is the subject of our current work in progress. Clearly, we could also have cho- 
sen another standard set of acceptance probabilities conforming to the detailed 
balance conditions, e.g. we could accept a transition 7^ 1— > 73+^^ := 6 with 
probability min ^1, exp(7i^''^'(75) — and a direct check of the de- 

tailed balance conditions, based on Proposition 1, would show that the law of 
A^l^'^^ is invariant with respect to such a dynamics. However, in this dynamics, 
in general we cannot efficiently control the size of the region affected in a single 
update. 

Versions of the disagreement loop birth and death dynamics can be easily 
constructed which leave invariant the distributions of the polygonal fields -^^I'g ' 

(•^Dl'biack' "^^Iwhite) '^i^h empty (black, white) boundary conditions respectively. 
To this end, we modify accordingly the dynamics (DL : birth[a, /3; a, b]) and 
(DL : death[a, /3;a, b]) by discarding all the updates which make the contour 
collection 7^ hit the boundary, and for the monochromatic black or white 
boundary condition, in addition, upon an update we do not pick the colour- 
ing by random but we choose the unique one compatible with the bound- 
ary condition. Denoting the so constructed dynamics by (DL0 : ...[a, /3; a, b]), 
(DLbiack : ...[a,/3;a,b]) and (DLwhite : ...[a,/3;a,b]) respectively, we immedi- 
ately conclude the following corollary from Theorem 1 

Corollary 1 For each a > 0, 6 > and a+a > 0, /3+b > 0, the law of the Gibbs- 

modified Arak process A^I^I^\A^I^'^^1^^^^, A^l^l^^^^^^) is the unique invariant distribu- 
tion distribution of the dynamics (DL0 : ...[a, /3; a, b]) /(DLbiack : ...[a. /?; a, b]) 
or (DLwhite '■ ...[q;, /3; a, b])) respectively]. The resulting stationary processes are 
reversible. For any initial distribution of 70 the laws of the random polygonal 
fields 7s converge in variational distance to the law o/^^jg' Mij'i'biack' -^Diwlite 
respectively], as s 00. 
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We believe that a very similar dynamics could be used to simulate length- 
and area-interacting modifications of more general consistent polygonal Markov 
fields admitting the dynamic representation as discussed in Arak & Surgailis 
(1989,1991) and Clifford, Arak & Surgailis (1993). The only change would 
be an appropriate redefinition of the operations A®[-;-] and A®[-;-], and the 
resulting disagreement field would no more be a single loop. 

2.2 Contour birth and death dynamics 

As already mentioned, unlike the previous one, the dynamics discussed in this 
subsection is constructed in a much narrower setting, restricted to colourless 
contoTir configurations which do not hit the boundary, and it is meant to leave 
invariant the distributions of -4^ 10'. Recall from the discussion following (10) 
that in this setting we can take D to be an arbitrary bounded open set in 
W^, with piecewise smooth boundary, we do not need convexity. The approach 
developed in this section leads to a simulation algorithm discussed in Section 
4 below, which, though perfect, seems to be practically infeasiblc due to non- 
constructive description of the intensity measure of contour births. However, 
its value lies in that its infinite volume extension provides important theoretical 
information about the thermodynamic limit yielding in particular the 

uniqueness of the thermodynamic limit for /? large enough. Observe that the 
dynamics constructed in this section could be in principle also used directly for 
Metropolis sampling, yet the previous disagreement loop dynamics seems much 
better suited for this particular purpose. 

To proceed, we consider the space Cd consisting of all closed polygonal con- 
tours in D which do not touch the boundary dD. For a given point configuration 
X := {xi, wc denote by Cd{x) the family of those polygonal contours in 

Cd which belong to r£)(a;), i.e. whose vertex set coincides with x. We construct 
the so-called free contour measure ©d on Cd by putting for C C Cd measur- 
able, say, with respect to the Borel cr-field generated by the Hausdorff distance 
topology. 



with the Hamiltonian $ as in (2), with Fin(£)) standing for the family of fi- 
nite point configurations in D and where v* is the measure on Fin(£>) given 
by di'*{x) := dxi...dxn- In order to provide an alternative line- rather than 




(11) 
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point-based expression for Qd, for a given finite configuration (l) := {h, ■■■,ln) 
of straight lines intersecting D denote by Cd{1) the family of those polygonal 
contours in Cd which belong to Td{1)- Then we have, see e.g. (3.8) in the proof 
of Theorem 1 in Arak, Clifford & Surgaihs (1993), 

Qd{C)= [ Y. exp(-21ength(^))d/i*((0) (12) 

JFin{L{D]) g^cnCnil) 

with Fin(L[i?]) standing for the for the family of finite line configurations inter- 
secting D and where /i* is the measure on F'm(L[D]) given by (i/i*((/i, In)) '■= 
dii(li)...dij{ln) with /X defined in the discussion preceding (1). 

For /3 s R we consider the exponential modification 6^' of the free measure 
©D, given by 

Q^§{de) :=exp{-l3length{e))QD{de). (13) 

It is easily seen that the total mass 6^'(Cd) is always finite. Indeed, using 
(12), taking into account that the length of a line segment in D can be at most 
diam(Z?) and recalling that, by standard integral geometry, M := /i({Z \ IDD ^ 
0}) < length(9conv(D)) we conclude that 

e^\CD) < f; ''Mm die^miP)) ^ 

fe=0 

exp[length(a conv(£))) exp(|/3| diam(£)))] < oo. (14) 

Let VaW] be the Poisson point process on Cd with intensity measure 0^'. It 
then follows directly by (11), by the point-based representation (3) and by (8) 
that for all /3 G K for which the partition function Eexp ^— ?l!^''''(^£)j0)^ in 
(10) is finite (in fact, wc show that this holds for all /3 G M in Corollary 2 below), 
the polygonal field -A^jg' coincides in distribution with the union of contours in 
■p„[3] conditioned on the event that they are disjoint so that 

( 



^K^)=^ u ^ 



(15) 



where the conditioning is well defined in view of (14). In particular, taking into 
account (1) and (12), we have for all /? where (10) makes sense 

P (Ve.e'ep^i^] 9 ^ 9' ^ 9 f^e' = ib\ ^ exp(-[2 + (3] length((5)) 
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with r£)|0 standing for the family of admissible polygonal configurations in D 
which do not touch dD. It easily follows that the law of A^^'l^ is invariant 
and reversible with respect to the following contour birth and death dynamics 
(7s)s>o on r£,|0. 

(C : birth[/3]) With intensity e^§{de)ds do 

• Choose a new contour 9, 

• If ^ n 7s = 0, accept 6 and set js+ds '■= 7s U 9, 

• Otherwise reject 9 and keep 7s+a!s := 7s, 

(C : death[/3]) With intensity 1 for each contour 9 € remove 9 from 7s 
setting 7s := 7s \ 9. 

It is worth noting that, should we accept all the new-coming contours without 
the disjointness test in the above dynamics, we would get the Poisson contour 
process V^m as the stationary state. 

Observing that the process 7s constructed above spends a non-null fraction 
of time in the state and using a standard coupling argument we are led to 

Theorem 2 The law of the Gibbs-modified Arak process .^^^^ is the unique 
invariant distribution of the dynamics (C : birth[/3]), (C : death[/3]). The re- 
sulting stationary process is reversible. For any initial distribution of 70 the 
laws of random polygonal fields 7s converge in variational distance to the law of 
A^jjI^' as s —> 00. 

All our results in this section were conditional on the partition function in (10) 
being finite. We claim here that this holds for all /? e K. Indeed, since (3 = 
clearly satisfies this condition as corresponding to the basic Arak process Ad\^, 
Theorem 2 can be used for (3 = 0. The dynamics (C) above implies that the 
empty-boundary Arak process AD\ii) is stochastically dominated (in the sense of 
inclusion) by the union of contours in T^Od, see Corollary 5 below. In particular, 
by (14), for all a,/? e M 

Eexp(-H^'''l(i^|0)) <exp(|a|A(£)))Eexp ^ length(^) 

exp{\a\AiD))cxpi-QDiCD))oxp {el^^HCo)) < 00. 
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By an appropriate redefinition of 0_d admitting edges chopped off by the bound- 
ary, the same argument can be repeated for Aoiti replaced by Ad- Thus, we 
have proven 

Corollary 2 For each bounded open domain D C both the partition func- 
tions Mexp (^—'HJ-^'^^Ad)^ in (6) and Eexp ^— W^''''(^£)|0)^ in (10) are finite 
for all a, /3 G K. 

3 Thermodynamic limit 

The purpose of this section is to define the notion of thermodynamic limit for the 

considered polygonal fields and to establish its existence (cf. Surgailis (1991)). 

For a smooth closed simple (non-intersecting) curve c in D by the trace of 
a polygonal configuration 7 on c, denoted in the sequel by 7 A c, we mean the 
knowledge of 

• intersection points and intersection directions of 7 with c, 

• colouring of points of c. 

This concept can be formalised in various compatible ways, yet we keep the 
above informal definition in hope that it does not lead to any ambiguities while 
allowing us to avoid unnecessary technicalities. For convenience we assume that 
no edge of 7 is tangent to c, which can be ensured with probability 1 in view of 
the smoothness of c. 

Fix Q!,/3 G R. In view of the Gibbsian representations (1), (3) and (6) we 
easily check that for each c as above and with 9 standing for a trace on c there 
exists a stochastic kernel .A|"j'^'(-|^) with the property that 

C,^,jA^n''\A^n''^^c = e)=C,^,^^^^ 

for all bounded open and convex D 13 Intc and for bd £ {black, white}, where 
^int c denotes the law of the argument random element restricted to Int c (the 
interior of c). 

Consider the family T^2 of whole-plane admissible polygonal configurations, 
determined by (PI), (P2) and (P4) ((P3) is meaningless in this context) 
and by the requirement of local finiteness (any bounded set is hit by at most a 
finite number of edges). Let be the corresponding collection of black-and- 
white coloured whole-plane admissible polygonal configurations. It is natural to 
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define the family Q^A^"^'^^) of infinite volume Gibbs measures (thermodynamic 
limits) for ^["■'^l as the collection of all probability measures on with the 
accordingly distributed random element A satisfying 

Antc(i|iAc = ^)=4^fJ(-|0). (17) 

In addition, wc shall consider the family GriA^"'^^) of isometry invariant mea- 
sures in Q{A^"'^^^). Using an appropriate relative compactness argument much 
along the same lines as in Schreiber (2004) we will readily get the existence of 
at least one isometry-invariant thermodynamic limit for each (3 > 0. 

Theorem 3 For all a gR and /? > 0, the family Gt{-^^"'^^) is non-empty. 

Note that for a = and /? large enough this statement follows also by Theorem 
in Surgailis (1991). 

In the sequel, we will establish certain uniqueness results for the thermody- 
namic limit in the low temperature region within a particular class of infinite- 
volume measures without infinite contours. However, we do conjecture that 
for a = outside this class there exists an infinite number of extreme infinite- 
volume phases breaking both the rotational and translational symmetries. We 
briefly and informally sketch their tentative construction. For the increasing 
sequence of squares (— n, n)^, n = 1,2, ... we consider a sequence of boundary 
conditions arising by requiring that a large number C{n) of edges hit the left- 
hand side of (— n, n)^ (with the intersection points located more or less uniformly 
over the edge), the same number of edges intersect the opposite right-hand side, 
but no edges hit the upper and lower sides. We believe that by choosing an 
appropriate growth rate for C(n) we can assure that the resulting sequence of 
polygonal fields on (— n, n)^ is uniformly tight (e.g. in the topology discussed 
in the proof of Theorem 3) and the accumulation points of this sequence are 
thermodynamic limits for ^1°''^] with infinite number of infinite left-to-right 
polygonal chains. Moreover, the expected number of such chains hitting a disk 
of radius 1 should exhibit untempered growth to infinity with the distance of 
the centre of the disk from the origin. We conjecture that such untempered 
thermodynamic limits should exist even for ,/? = where, in the language of the 
dynamic time-space construction of the basic Arak process, one could, roughly 
speaking, have an infinite-density cloud of particles born at the time -co. Such 
constructions are possible due to the fact that, under very rapid edge density 
growth with the distance from the origin, one can enforce the situation where 
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the influence of the boundary conditions on d{—n, n)^ competes on equal rights 
or even dominates the stabilising bulk effects within (— n,n)^. Clearly, such 
phenomena cannot show up in the stationary regime, see Schreiber (2004) for a 
discussion. 

4 Perfect simulation from thermodynamic limit 
and exponential mixing 

The purpose of the section is to study the contour birth and death dynamics 
of Subsection 2.2 in context of the perfect infinite- volume simulation scheme as 
developed by Fernandez, Ferrari & Garcia (1998,2002). This approach is valid 
only for sufficiently large (3. It yields a perfect algorithm for simulating ther- 
modynamic limits in finite windows and it allows us as well to conclude certain 
uniqueness and mixing results for the thermodynamic limit in low temperature 
regime. 

To this end, we observe first that for all bounded open sets D with piecewise 
smooth boundary the free contour measures as defined in (11) arise as the re- 
spective restrictions to Co of the same measure 9 = 9r2 on C := U^i C(-n,ra)2, 
in the sequel referred to as the infinite volume free contour measure. Indeed, 
this follows easily by the observation that restricted to coincides with 
for D2 C Di. In the same way we construct the infinite- volume exponen- 
tially modified measures Ol'^' = ojjfj- The following result, which is related to 
the Lemma in the Appendix of NichoUs (2001), will be crucial for our further 
purposes as stating exponential decay of the measure Q^^^ with respect to the 
contour size. 

Lemma 1 For (3>2 we have 

e^l^\{e I dx e Vertices(6'), length(6') > R}) < 47rexp(-[/3 - 2]R)dx. 

(18) 

Moreover, there exists a constant e > such that, for /3 > 2, 

el''] ({6* I € Int6i, length(6') > R}) < exp(-[/3 -2 + e/2]R + o{R)). 

(19) 

We note that, in view of (15) in Section 2.2, a standard Peierls-type argument 
can be applied to conclude from Lemma 1 that there is no infinite contour 
nesting for ^['^'''1 whenever (3 >2. 
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The approach of Fernandez, Ferrari & Garcia (1998, 2002) speciahsed for our 
purposes rehes on the following graphical construction, briefly sketched below, 
see ibidem for further details. Choose /3 > 2 large enough, as specified below. 
Define J-{C) to be the space of countable and locally finite collections of contours 
from C, with the local finitcness requirement meaning that at most a finite 
number of contours can hit a bounded subset of M? . Observe that J^{C) C 
(there is no equality since J-{C) contains only bounded closed contours while 
riK2 also admits infinite polygonal chains). On the s-time-space M x T{C) we 
construct the stationary unconstrained (free) contour birth and death process 
(^3s)sgR with the birth intensity measure given by O^'^l and with the death 
intensity 1. Note that unconstrained or free means here that every new-born 
contour is accepted regardless of whether it hits the union of already existing 
contours or not, moreover we admit negative time here, letting s range through 
R rather than R+ . Observe also that we need the birth measure to be finite 
on the sets G C | 6'n A 0} for all bounded Borcl A C in order to have the 
process (es)sgR well defined on M x J^(C). By Lemma 1 this is ensured whenever 
/3 > 2. It is easily seen that, for each s G M, £>s coincides in distribution with 
the whole-plane Poisson contour process Vq[i3\ . 

To proceed, for the free process {Qs)se^ ^"^c perform the following trimming 
procedure. We place a directed connection from each s-timc-spacc instance of 
a contour showing up in (gs)sgR and denoted by x [so,si), with stand- 
ing for the contour and [so,si) for its lifespan, to all s-time-space contour in- 
stances 6' X [sq,s'i) with 6' r\ 9 ^ %,s'q < sq and s'^ > sq- In other words, 
we connect 6 x [sq, si) to those contour instances which may have affected the 
acceptance status of ^? x [,so,.si) in the constrained contour birth and death dy- 
namics (C) as dicussed in Subsection 2.2. These connections yield directed 
chains of s-time-space contour instances, we call them the ancestor chains in 
the sequel. Following Fernandez, Ferrari & Garcia (2002) the union of all an- 
cestor chains stemming from a given contoiir instance is referred to as its clan 
of ancestors. Using Lemma 1 combined with a general technique of stochastic 
domination by subcritical multitype branching processes as discussed in detail 
in Fernandez, Ferrari & Garcia (1998,2002), for (3 large enough we can ensure 
that all such clans of ancestors are a.s. finite and that a single clan size has 
exponentially decaying tail [i.e. the probability that the clan size exceeds R 
is of order 0(exp(— ci?)) for some c > 0]. In this case we can uniquely deter- 
mine the acceptance status of all the clan members: contour instances with no 
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ancestors are a.s. accepted, which automatically and uniquely determines the 
acceptance status of all the remaining members of the clan by recursive appli- 
cation of the inter-contour exclusion rule. Discarding the unaccepted contour 

instances leaves us with an s-time-space representation of a stationary evolu- 
tion on J-{C) C r]g2. The graphical construction and the argument in 
Fernandez, Ferrari & Garcia (1998, 2002) specialised to our setting yield 

Theorem 4 Choose (3 > 2 large enough so that all the ancestor clans in the 
above graphical construction are a.s. finite and a single clan size exhibits expo- 
nentially decaying tail. Then 

1. the J-" (C) -valued process (7s)s>o given above is well-defined, stationary and 
reversible, 

2. the stationary distribution £(70) on T{C) C is isometry invariant and 
belongs to GriA^^-'^^), 

3. the dynamics of {js)ses. is an infinite-volume extension of the contour 
birth and death dynamics (C) as introduced in Section 2.2, i.e. {'ys)seM is 
a Markov process on !F{C) with the infinitesimal generator 

[mF]{r^) := / [F{n U {6}) - F{n)]dQ^^\e) + Y.[F{n \ {0}) - F{v)] 
•^c eer, (20) 

for T] G J^{C) and bounded F : J-{C) M such that Firf) depends only on 
rjD D for some bounded convex set D, 

4- (7s)seR exhibits exponential s-time-space a-mixing in that there exists c > 
such that 

sup |P(£:i n £2) - P(fi)P(f2)| < e-'=['^'"*(^'^)+'*'"*«''0'*il't''o'''il)l 

£ie3B(n=,l)x[so,oil 

whenever dist{x,y) is sufficiently large, with '<sb(x,i)x[sq,si] standing for 

the (j-field generated by the restriction o/(7s)sgR to the s-space-time region 
B(x, 1) X [so, si], where B{x, 1) is the disk of radius 1 centred at x £ M^, 

5. consequently, the stationary distribution £(70) exhibits exponential spatial 
a-mixing. 
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It is worth noting that even if (3 is not large enough to ensure a.s. finiteness 
of ancestor clans, a weaker version of the above graphical construction can be 
provided as soon as the birth intensity measure el^l is finite on {6 & C \ ef\A 
0} for all bounded A C.M? , which is the case whenever /3 > 2 by Lemma 1. To 
this end we restrict the s-timc to M_|- and choose an initial condition, which is an 
J^(C)-valued random element independent of the free birth and death process of 
the graphical construction. The birth and death process here is also restricted to 
positive times in that there are no contours born or alive before the s-time 0, in 
other words the birth and death process starts with the initial state at s-time 0, 
consequently it is no more stationary. In this context the local finiteness of G^^l 
allows us to conclude that for each contour instance x [sq, si), so, si > 0, the 
expected cardinality of its ancestor clan extending down to the s-time is finite, 
consequently the clan is a.s. finite (note that it could extend to an infinite clan 
through negative s-times in the original graphical construction). Thus, with the 
initial state given, the acceptance status of each contour instance is uniquely 
determined by the intcr-contour exclusion rule. This leads us to 

Corollary 3 With /3 > 2, for each T{C)-valued initial condition 70 there exists 
a Markov process (7s)s>o on T{C) with infinitesimal generator given by (20). 

In the remaining part of the present section wc will not use Corollary 3 and, 
unless otherwise stated, we shall assume that (5 stays within the region of valid- 
ity of the original graphical construction preceding Theorem 4. We denote by 
/^I^l the infinite-volume stationary distribution £(70) arising in this graphical 
construction. Observe that the fact that /i^^l is concentrated on J-{C) means 
that it contains no infinite polygonal chains - all the contours are bounded and 
closed. Below, we show that, with the assumptions of Theorem 4, jj^^ is in 
fact the unique element of Q{J^^'^^) concentrated on !F{C), although we con- 
jecture that Q{J^^'^^) is non-empty as argued in Section 3. To proceed with 
our argument we consider finite-volume versions of the above graphical con- 
struction, with the infinite- volume birth intensity measure G^'^' replaced by its 
finite volume restrictions 0^' for bounded and open D with piecewise smooth 
boundary. Clearly, the graphical construction yields then a version of the finite- 
volume contour birth and death dynamics (C) as discussed in Subsection 2.2. 
For each D denote the resulting finite- volume stationary process on T{Cd) by 
{l?)sew.- Write also {g'^) for the corresponding free contour birth and death 
process. Note that this finite-volume construction is valid for all /3 € R, even 
though in this section it is only used for (3 as in Theorem 4. In view of Theorem 
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2 we see that coincides in distribution with ^^jg for all s € M. Moreover, it 
is easily seen that coincides in distribution with V^w for all s € K. Prom the 
construction, Lemma 1 and the general theory developed in Fernandez, Ferrari 
& Garcia (1998, 2002) it follows that 

Proposition 2 With (3 as in Theorem 4 the finite-volume graphical construc- 
tions for different Z) C and the infinite-volume graphical construction can be 
coupled on a common probability space so that there exists c > with 

P (7f 1 n B{x, 1) 7^ 7f " n B{x, 1)) < exp(-cmin(dist(a;, dDi), dist(a;, 5£>2))) 

for bounded Di,D2 C R^, for x sufficiently far from dDi and dD2 and for all 
s G M. Moreover, 

P (7f n B{x, l)^j^n B{x, 1)) < exp(-cdist(a;,5L')) 

for bounded D CM?, for x far enough from dD and for all s G R. 

Taking into account that, by the construction and by the results of Section 2.2, 
7^ coincides in distribution with j^'^^^ , and that for each contour collection in 
Tic) every bounded region can be surrounded by a smooth curve which does 
not hit any of the contours, we can use the Markov property of the considered 
polygonal fields combined with Proposition 2 to conclude the claimed property 

Corollary 4 For (3 as in Theorem 4 the measure is the only element of 
g{A^°'f^^) concentrated on T{C). 

For (3 as in Theorem 4, using Lemma 1 we easily conclude that the number of 
contours in qq surrounding a given point is a.s. finite. Consequently, the number 
of contours surrounding a given point in 70 is a.s. finite as well, whence there is 
no infititc contour nesting. Thus, wc observe a unique infinite connected region 
surrounding finitely nested contour collections. Colouring this region black or 
white gives rise to two distinct phases, respectively black- and white-dominated. 
There are no other extreme phases without infinite chains in the coloured model, 
because their corresponding colourless contour ensembles have to coincide with 

The last important conclusion of the graphical construction, based on the 
above-made observations that almost surely 7s C gs,j^ C and that 7^ 
^[o,/3]^^r> A^°'l^^\gs Vsi0] and £»f r-.^ -p^^^^^ the following stochastic 
domination statement 
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Corollary 5 The Poisson contour process Vqi/i] stochastically dominates (in 
the sense of inclusion of contour collections) the polygonal field A^^'^^. Likewise, 
for each bounded D with piecewise smooth boundary, the Poisson process V^m 

stochastically dominates the finite-volume polygonal field A^^'^^ . 

5 Proofs 

5.1 Proof of Theorem 1 

In order to provide a geometrical and intuitive proof of the theorem we construct 
an auxiliary model. For r > define A^^'^'^^ to be the Gibbsian modification of 
Ad with the Hamiltonian 

:= r-'pA ([7 +m B{r)] n D) + aA((black[7] +m B{r)) n D), 

(21) 

with +M standing for the usual Minkowski addition and with B{r) denoting 
the radius r disk in M?, centred in 0. It is easily seen that, for each 7 e Fn, 

lunH^S'''''\^)=n^S'%) (22) 

so that H^'''-''' is an approximation of H^'^' for small r. Take nl^+'^l , n[''"'(''+^)] , 
n["l and nl'" ^''1 to be independent homogeneous Poisson point processes on 
D, jointly independent of Ad, with respective intensities a + a,r~^{l3 + b),a 
and r~^b. We claim that A^^''^'^^ coincides in distribution with Ad conditioned 
jointly with n["+'^l,n[''"'(''+'')l, nl"^] and nl''"'*'] on the event £[a, P;a,b;r] that 
the following conditions are simultaneously satisfied 

. U^--'(P+b)]n[-f+MB{r)]=$, 

• nl'^""] C [7+MB(r)], 

• n[«+«I n [black(7) +M B{r)] = 0, 

• nW C [black(7) +M -B(r)] = 0, 
so that 

£(^^■^^'■1) = C (AD\£[a,P;a,b;r]) . (23) 
Indeed, for a given 7 G the probability of the event £[a, (3; a, b; r] is 
F{£[a, (3; a, b; r] I7) = exp (-r'^ [/3 + b]A{['y +m B{r)] n D)) 
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exp {-r~^b[A{D) - ^([7 +m B(r)] n 7^)]) 
exp(-[Q; + a]A([black(7) +m B{r)] n D)) 
exp {~a[A{D) - A([black(7) +m B(r)] n £))]) = 
exp (^-H^^'^-''\^)) exp (-[a + r-'b]A{D)) , 

which yields (23) by the definition of A^^'^'^^ . 

To proceed, we construct an auxihary Markovian dynamics which leaves 
invariant the joint distribution of of ic, n[''"'('^^+*)l , n["+"l , hI''"'''] and H^, 
and makes the resulting stationary process reversible. To this end, set 

70 := Ad, ~ := < := H^, ^ := nl*-] 

and let the quintuple (7-5, tt", 7rf , 7r°, 7rg)s>o evolve according to the following 
rules, applied independently for each component, 

(Auxl) 7s evolves according to (DL : birth) and (DL : death), 

(Aux2) 

■""^ ) ■"'f ! ""s TTg evolve according to a birth and death process with 
death intensity 1 and with birth intensities a + a, r^^(/3 + b), a and r~^b 
respectively. 

The above invariance and reversibility statements follow as direct consequences 
of Proposition 1. Thus, we conclude that the joint distribution of {Ad,^^^ 
jj[a+a]^ jj[r ^fej^nl"]) Conditioned on the event f [a, /3; a, b;r\, is invariant and re- 
versible with respect to the following Markovian dynamics, arising from (Auxl) 
and (Aux2) by adding an appropriate acceptance test to be passed only by ad- 
missible updates: 

(Bl) Choose an update (S, ^ , f?^ e^.^e'^) for (7,+^., 7rf+,„ tt^^^J 

according to the rules (Auxl), (Aux2). 

(B2) Accept the update, setting 

(7.+d.,<+d„<d..<+d..'^'+dJ ■■= {s,6'',e^6\e''), 

provided the following conditions are satisfied 
• 0" C [6+MB{r)], 
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• 61" n [black((5) +M B{r)] = 0, 
. 61" C [black((5) +MS(r)], 

(B3) Otherwise discard the update, keeping 

Consequently, in view of (23), the first component 7,, under the above stationary 
dynamics (Bl — 3), with the initial distribution at s = given by the joint law 
of (ir,,n[''"'(^+*)l, n[«+»l, nl'-"'*'!, nW) conditioned ontheevent£[a,/3;a,6;r], 
coincides in distribution with A^^'^'^^ for all s S Moreover, the conditional 
distributions of the remaining components given 7^ are also readily determined. 
Indeed, tt" is just a homogeneous Poisson point process on D \ [black(7s) +m 
B{r)] with intensity a + a, 7rf is an intensity r^^(/3 + 6) homogeneous Poisson 
point process on Z) \ [7 +m B{r)], tt" is a homogeneous Poisson point process 
on black(7s) +m B{r) of intensity a while tt^ is a homogeneous Poisson point 
process on 7^ +m B{r) with intensity r~^b. All four components tt" , 7rf , 7r° , tt^ 
are jointly independent given 7^. Consequently, we observe that if we integrate 
out the Poisson components and tt'', the polygonal field component 

7s turns out to evolve according to the following dynamics (see Subsection 2.1 
for the notation): 

(DL : birth[a, /3; a, b; r]) With intensity [wdx + n{dx)\ds do 

• put (5 := 7s © X, 

• construct 5 by randomly choosing, with probability 1/2, either of the 
two possible colourings for 5, 

• accept 5 with probability 

exp (-[a + a\A ([black(5) +m B{r)\ \ [black(7,) +m S(r)])) 

exp (-r-M/3 + h]A{[5 +m B{r)] \ [7^ +m B{r)])) 
exp (^-aA ([black(7s) +m B{r)] \ [black((5) +m 5(r)])) 

exp (-r-i6^([7, +m B{r)] \ [5 +m B{r)])) = 
exp [-aA ([black(5) +m B{r)] \ [black(7s) +m B(r)])) 
exp {-I3r-^A{[5 +m B{r)] \ [7^ +m -B(r)])) 
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exp 



(-aA ([black(,5) +M -B(r)] A[black(7,) +m 
exp {-br-^A{[6 +m B{r)]A[js +m B{r)])) , 



• if accepted, set ^s+ds '■= 



S, otherwise keep js+ds '■= 7; 



(DL : death[a, /3; a, b; r]) For each birth site x in 7^ with intensity 1 do 
• put (5 := 7s X, 



Thus, the distribution of A^^'^'^^ is invariant and reversible with respect to the 
above dynamics. Moreover, it is easily seen that the acceptance probabihties 
in the rules (DL : birth[Q:, /3; a, b; r]) and (DL : death[Q!, /?; a, b; r]) converge 
to these in (DL : birth [a, /3; a, b]) and (DL : death[a, /3; a, b]) as r ^ 0. Tak- 
ing into account (22) and letting r ^ we get by a standard continuity ar- 
gument that A^^'^^ is invariant and reversible with respect to the dynamics 
(DL : birth[a,^;a,b]) and (DL : death[a, /?; a, b]). 

To complete the proof of Theorem 1 it suffices now to establish the remaining 
uniqueness and convergence statements. These follow, however, along the same 
lines as in Proposition 1, by the observation that, in finite volume, regardless of 
the initial state, the process 7^ spends a non-null fraction of time in the state 
'black' (no contours, the whole domain D coloured black) and by a standard 
application of coupling argument. The proof is complete. □ 



• construct d by randomly choosing, with probability 1/2, either of the 
two possible colourings for 6, 

• accept d with probability 




• if accepted, set ^s+ds '■= 



6, otherwise keep 'js+ds '■= 7< 
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5.2 Proof of Theorem 3 



Following the ideas of Schrcibcr (2004) it is convenient to consider the family 
of admissible configurations in the plane embedded into the space G]s2 of 
locally finite non- negative Borel measures on M?, by identifying a configuration 
7 e with the measure 

M^{U) := length(7 nU) + yl(black(7) nU) + N{j n U) (24) 

for Borel [/ C R^, with N{'jr)U) standing for the number of vertices of 7 falling 
into U. Endow the space Gr2 with the vague topology defined as the weakest one 
to make continuous the mappings H J fdfj, for all continuous / with bounded 
support. Observe that in general <^ T]g,2 for G C due to the presence of 
edges chopped off by the boundary. Therefore, in order to have our embedding 
defined also for finite- volume configurations, we agree to put M^(D^) :— for 
all 7 e r£). Note that only the internal vertices of finite- volume configurations 
are counted in N{-). 

Consider the sequence ((— n, n)^)^i of growing open squares in R^. By the 
properties of the basic Arak process, sec Section 4 in Arak & Surgailis (1989) 
and Section 2.1 in Schreiber (2004), it immediately follows that there exists a 
finite constant C with 

EM .-[0,01 {{-n,nf) <CA{{-n,nf) (25) 

■^(-n,„)2 

for all n > 1. We will show that the above conclusion can be extended for 
arbitrary a G M and /3 > in that there exists C'"'^! < 00 with 

EM,-[.,si ((-n,n)2) < (7["''^U((-n,n)2). (26) 

Below, we assume without loss of generality that a > 0, which can be done in 
view of the colour-flip symmetry. Observe first that, in view of (6), 

with H^^^^y as in (7). Consequently, taking into account that the area term 
in the Hamiltonian 'H^(^^^^2 is bounded by aA{{—n, r?.)^) and that the Hamilto- 
nian is always positive, we conclude by (25) that also the expectation of the edge 
length term in the Hamiltonian admits an area-order upper bound. It remains 
to show that this is also the case for the number of vertices - we sketch the argu- 
ment omitting standard technical details. To this end, we take advantage of the 
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dynamic representation (as discussed in the introduction of this paper and in 
Section 4 of Arak & Surgailis (1989)) to conclude that for the basic Arak process 
•^(-n,Ti)2 the number of internal left-extreme vertices (with the corresponding 
sharp angle lying to the right of the vertex) is Po(7r^((— n, n)^)), where Po(t) 
stands for Poisson-distributed random variable with mean t. The same applies 
for the number of internal right-extreme, upper-extreme and lower-extreme ver- 
tices (recall that we do not count the boundary vertices here). Consequently, 
the overall number of internal vertices N{A(-n,n)'^) is stochastically bounded by 
4Po(47TO^) and has its mean of area order, not greater than 167rn^. In view of 
the representation (6) and taking into account that the Hamiltonian '^("'f' „)2 
is always positive since a > 0, we conclude that, for all iiT > 0, 

^ ^""'"^ ^ " Eexp(-HKl,.(iKl ,.)) 

Recall that Poisson distributions exhibit supcrcxponentially decaying tails 

P(Po(47rn2) >K)< exp (^-^ log (^-^^^^ ■ ^ > 647rn^ 

see Shorack & Wellner (1986), p. 485. Moreover, the negative logarithm of the 
denominator in (27) exhibits at most area-order growth, which is due to the 
easily verified finiteness of the free energy density for T^I"''^] 

Consequently, the required area-order bound for E7V(^'"'^'^^^2) follows now from 
(27) by a direct calculation. This completes the verification of (26). 

To proceed with the proof of the theorem, consider the sequence (M^""^')^^ 
of (jR2 -valued random elements with laws given by 

mt^^) ■■= j4r^ [ [ C ([t, o R^jM^i..,, ] dxdct>+ 
47r(2n)^ i[0,2,r) ■/(-n,n)2 V -^(-n.nj^y 



^^—2 / / £([EoT^oi?^]M,[c.^i ]dxd(j), 



(28) 



where Tx stands for the standard translation operator TxH{U) := ^{U + x) while 
ii,^, G [0, 27r) is the rotation by angle (f) around and S is the reflection with 
respect to some fixed axis passing through the origin. By (26) it follows that 

< 00 
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for all bounded U C IR.^. Applying Corollary A2.6.V in Daley & Vcrc- Jones 
(1988) we conclude that the sequence of random measures {M^'^^)'^^-^ is uni- 
formly tight in Gr2 and, consequently, it contains a subsequence converging in 
law to some M^c corresponding to a whole-plane polygonal field a!S>'^^ ■ In view 
of (28) it is clear that 

which completes the proof of the theorem. □ 
5.3 Proof of Lemma 1 

By the definition (13) of the /?-tilted contour measure ©[''I it is enough to es- 
tablish the assertion of the lemma for the henceforth assumed case /? := 2. In 
order to establish (18) define the continuous-time random walk Zt in M? with 
the following transition mechanism 

• between critical events specified below move in a constant direction with 
speed 1, 

• with intensity given by 4 times the covered length element update the 
movement direction, choosing the angle (f) E {0, 2tt) between the old and 
new direction according to the density | sin(^)|/4. 

We start the random walk Zt at a given point x and with a given initial velocity 
vector. Moreover, we choose the loop-closing angle (j)* G (0, 27r) according to the 
density | sin(0)|/4 and we draw an infinite loop-closing half-line I* starting at x 
and forming the angle (jf with the initial velocity vector. Let Zt be the random 
walk Zt killed whenever hitting its past trajectory or the loop-closing line I*. The 
directed nature of the random walk trajectories as constructed above requires 
considering for each contour 9 two oriented instances 9^ (clockwise) and 9^ 
(anti-clockwise). We claim that for x e and 9 & C with x € Vertices(6) we 
have 

87rda;e"'"'"'sth(e-) p reaches 1* and the resulting contour falls into d6'~*) 

= ePl(d0), (29) 

where e* stands for the last segment of 9^ counting from x as the initial vertex, 
which is to coincide with the segment of the loop-closing line 1* joining its 
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intersection point with Zt to x. Clearly, the same relation holds then for , 
hence adding versions of (29) for and , which amounts to taking into 
account two possible directions in which the random walk can move along 
will yield 261^1(^61) on the RHS. The relation (18) will easily follow by using 
the trivial upper bound 1 for the probability on the LHS of (29). 
To establish (29), we observe that the probability element 

P {^t reaches I* and the resulting contour falls into dB~^^ 

is exactly 

k 

%XM]({(M*)Mn^gd.}) -P(-41ength(^ \ e*)) [J ^V.W ^^^^ 

where ei, . . . , are all segments of d including e*, while l\ej\ stands for the 
straight line determined by e,;. Indeed, 

• the prefactor [4[/i x /;i]({(?,r) ^ cZa;})]"-^ comes from the choice of 
the lines containing respectively the initial segment of [counting from 
a;] and T, as well as from the choice between two equiprobable directions 
on each of these lines, 

• for the remaining segments we use the fact that, for any given straight line 
Zo, I i n ^0 S (i^, Z(Z, Zo) € '^'/>}) = I sm.(^\dld4> with dl standing for the 

length element on ?o and with Z{lo,l) denoting the angle between / and 
/o, see Proposition 3.1 in Arak & Surgailis (1989) as well as the argument 
justifying the dynamic representation in Section 4 ibidem. Note that the 
direction update intensity was set to 4 to coincide with J^^ \ sm(f)\ = 4. 

To get the required relation (29) it is now enough to use (30), recall the definition 
of 6 and observe that [n x /z] {{{I, I*) \ iDl* e dx}) = 2iTdx as follows by standard 
integral geometry. This completes the proof of (18). 

To proceed, let Zf be the random walk Zt killed whenever hitting its past 
trajectory, but not when hitting the loop-closing half-line /*. Define 

^■■=-^^^J^^0S^if>T), (31) 

where f is the lifetime of Zt or, in other words, the first moment when Zt hits 
its past trajectory. The existence of the limit in (31) follows by a standard 
superadditivity argument, see Section 1.2 in Madras & Slade (1993), and in fact 
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£ can be regarded as the connective constant for the self-avoiding version of the 
random walk Zt, see ibidem. It is easily checked that e > since during each 
unit time of its evolution the walk Zt has a certain positive probability of hitting 
its past trajectory, uniformly bounded away from through time. To establish 
(19) observe that, as in the argument above, 

ePl({6» I dx G Vertices(6i), length(6») > i?}) < 

47rrfa;P [Zt survives up to time i?/2^ < 47rda;P(f > i?/2). (32) 
The required relation (19) follows now by (32), (31) and by the observation that 
e[2l({6' I € Int6i, length(6') > K\) < 

OO 

^ e[2l({6' I Verticcs(6') n [B(0, /c + 1) \ B(0, fc)] / 0, lcngth(6') > max(i?, fc)}). 
fe=o 

The proof is complete. □ 
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